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A  METHOD  OF  REDUCING  BEARING  ERRORS  IN  HORIZONTAL  LINEAR 
HYDROPHONE  ARRAYS  IN  SHALLOW  WATER 

by 

Richard  Klemm 


ABSTRACT 


Bearing  estimation  by  means  of  a  horizontal  acoustic  line  array  is  degraded 
seriously  by  the  multiple  effect  caused  by  dispersion  of  the  field  in 
normal  modes.  Error-free  bearing  estimates  can  be  achieved  only  by  means 
of  generalized  power  estimators  that  take  into  account  the  spatial 
statistics  (i.e.  the  covariance  matrix)  as  a  priori  information  about  the 
medium.  This  a  priori  information  is  unknown  in  reality.  It  has  to  be 
substituted  either  by  running  a  modelling  program  (e.g.  normal  mode)  or  by 
measurements  using  a  test  source.  Measurements,  however,  have  to  be 
carried  out^for  all  bearings  of  interest.  A  method  is  suggested  that  uses 
only  one  test  source  at  one  known  direction.  The  estimated  test  covariance 
function  is  used  to  calculate  a  beamformer  matrix  that  can  be  steered  in 
bearing.  Theoretical  investigations  show  that  significant  reduction  of  the 
bearing  errors  can  be  achieved. 
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envelope  of  approximate  spatial  covariance  function 
envelope  of  spatial  covariance  of  test  signal 
correction  exponent 
modal  amplitudes 

approximated  maximum  entropy  method 

bearing  angle 

bearing  of  target 

corrected  bearing 

look  direction  of  steering  matrix 

bearing  of  test  source 

correlation  decay  factor 

conventional  beamformer 

spacing  between  two  receivers 

bearing  error 

generalized  beamformer 

generalized  maximum  likelihood  estimator 

generalized  maximum  entropy  method 

steering  matrix 

steering  vector 

modal  wavenumbers 

wavelength 

number  of  modes 

maximum  entropy  method 

maximum  likelihood  estimator 

noise  covariance  matrix 

optimum  maximum  entropy  method 
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processor  output  power 

approximated  vector  of  phase  coefficients 

pressure 

vector  of  spatial  test-signal  phase  coefficients 
modal  phases 

observed  covariance  matrix 
range 

spatial  covariance  coefficient 

phase  correlation  between  n*'*1  and  m*'*1  mode 

covariance  coefficient  due  to  pilot  source  test  signal 

spatial  covariance  matrix  of  target 

target  signal 

trace  of  a  square  matrix 

width  of  equivalent  bearing  interval 

mid  of  equivalent  bearing  interval 
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INTRODUCTION 

This  paper  presents  a  method  of  reducing  the  bearing  errors  of  horizontal 
line  arrays  that  may  be  caused  by  the  multipath  effects  of  sound  propaga¬ 
ting  in  shallow  water.  In  a  waveguide  like  the  shallow-water  channel, 
sound  propagates  in  eigenfunctions  of  the  channel  (normal  modes),  each  of 
them  corresponding  to  a  wavefront  that  is  nonplanar  in  the  vertical  and 
arrives  at  a  certain  vertical  angle  at  the  receiver.  Even  if  different 
modes  are  completely  coherent  among  each  other,  the  received  wavefront  is 
in  general  nonplanar  also  in  the  horizontal.  If,  in  addition,  there  are 
phase  fluctuations  among  modes  due  to  some  randomness  of  the  channel  the 
observed  wavefront  is  random  and  nonplanar.  In  this  case  the  signal  has  to 
be  described  by  its  spatial  covariance  matrix  rather  than  by  a  vector  of 
spatial  field  samples. 

A  recent  theoretical  study  [lj  has  shown  that  such  well-known  resolution 
methods  as  conventional  beamforming,  maximum  likelihood  estimator,  and 
maximum  entropy  method  yield  serious  bearing  errors  for  all  angles  other 
than  broadside.  In  particular,  at  endfire  direction  the  bearing  offset 
may  amount  to  more  than  10°.  In  many  cases  more  than  one  peak  is  achieved 
as  array  response  to  a  point  source.  The  bearing  offset  and  the  number  of 
peaks  achieved  depend  on  the  actual  channel  and  source  parameters,  such  as 
the  water  depth,  frequency,  and  sound-speed  profile.  Furthermore,  they 
depend  on  the  array  aperture  and  the  resolution  method  used,  and  on  the 
positions  of  source  and  array. 

It  was  also  shown  in  [i]  that  these  problems  can  basically  be  overcome  by 
use  of  so-called  generalized  resolution  methods  that  make  use  of  the 
a  priori  knowledge  of  the  spatial  statistics  of  the  signal.  These  are 
generalized  in  the  sense  that  the  steering  vector  contained  in  most  of  the 
known  resolution  methods  is  replaced  by  a  steering  matrix,  which  is  a  set 
of  spatial  covariance  matrices  related  to  any  expected  signal.  Let 
R  =  S(60)  +  N  be  the  observed  spatial  covariance  matrix  consisting  of  a 
signal  and  a  noise  component.  Then  the  power  output  of  the  conventional 
beamformer  (CBF)  is,  in  the  mean, 

P((3)  =  h*(B)  R  h(6),  (Eq.  1) 


where  £  is  the  bearing  and  hjs)  is  the  steering  vector  with  the  elements 
h.j  =  exp(j2ird/X*i»cosB) ,  d  being  the  spacing  between  sensors.  Equation  1  is 
optimum  in  the  signal -to-noise  ratio  sense  if  the  signal  component  of  R  is 
deterministic,  i.e.  S  =  £(B0)  £*{ep),  and  if  h/eo)  =  s(B0),  which  means 
that  the  steering  vector  is  perfectly  matched  to'the  signal  wavefront.  If 
s( Bq )  is  random  the  steering  vector  may  be  replaced  by  a  steering  matrix 
H(B]  so  that  H(80)  -•  S(80).  Hence  the  generalized  beamformer  (6BF)  output 
power  becomes 


P(6)  =  tr(R  H(B)). 


(Eq.  2) 


Similarly,  the  maximum  likelihood  estimator  (MLE)  [2] 
P(3)  =  h*(B)  R'1  h(6)f 1 


(Eq.  3) 
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and  the  maximum  entropy  method  (MEM)  [3] 

P(8)  =  P0(e*R‘lh(B)h*(8)R'1e)'1,  (Eq.  4) 

with  e*  =  (1,0 . 0)  -  jo,  can  be  generalized  to  random  wavefronts 

P(8)  =  tr(R-1H(8) )-1  (GMLE) ,  (Eq.  5) 

P(8)  =  P0 ( e*R" 1 H ( 0 ) R" 1 e ) " 1  ( GMEM ) .  (Eq.  6) 


Equations  2,  5  and  6  represent  the  optimum  resolution  methods  if 
H(60)  =  S(80).  We  are,  however,  still  faced  with  the  problem  that  the 
steering  matrix  H(8)  is  unknown.  One  way  of  achieving  a  substitute  for 
H(b)  is  to  run  a  normal-mode  program  and  calculate  the  covariance  matrix  of 
expected  signals,  as  will  be  described  in  the  next  chapter.  However,  in 
this  case  the  performance  of  resolution  methods  depends  seriously  on  how 
well  the  predicted  signal  covariance  fits  the  observed  signal.  In  other 
words,  the  resolution  depends  on  how  far  it  is  possible  to  model  the  actual 
environment.  In  particular,  the  bottom  is  usually  unknown.  Nevertheless, 
this  method  has  to  be  investigated  experimentally. 

The  method  suggested  in  the  following  uses  a  pilot  source  with  known 
direction  to  transmit  a  test  signal.  This  is  used  to  estimate  the  spatial 
response  of  the  channel  to  a  point  source,  i.e.  the  matrix  H(Bt),  where  8t 
is  the  direction  of  the  test  source.  The  problem  is  investigated  theo¬ 
retically  by  using  a  normal -mode  program  [4]  to  model  the  channel  responses 
to  the  pilot  source  and  to  the  target  to  be  located. 
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1  THE  SIGNAL  MODEL 

The  SNAP  normal -mode  shallow-water  sound-propagation  model  [4]  gives  a  far 
field  solution  of  the  wave  equation  in  terms  of  the  eigenfunctions  of  the 
acoustic  waveguide 

M 

p(r,t)  =  exp(-joit)  E  a  exp(jk  r),  (Eq.  7) 

n=l  n  n 

where  p  is  the  pressure,  kn  are  the  modal  wavenumbers,  and  r  is  range. 

The  an  are  the  modal  amplitudes,  which  depend  on  the  positions  of  source 
and  receiver  and  several  channel  properties,  such  as  the  sound-speed 
profile  and  modal  attenuations.  The  values  an  and  kn  are  the  output  of 
the  SNAP  model  [4]  used  in  the  following.  The  expressions  kpr  determine 
the  modal  phases  between  source  and  receiver.  In  order  to  model  some 
randomness  of  the  boundaries,  an  additional  modal  phase  $n  is  introduced: 

M 

p(r,t)  =  exp(-juit)  E  a  exp[j(k  r)  +  <(>  )].  (Eq.  8) 

n=l  n  n  n 

This  may  be  zero  for  entirely  coherent  modes  or  more  or  less  random  when 
modelling  modal  phase  fluctuations.  Now  the  spatial  correlation  between 
two  points  of  the  field  r  and  r  +  d  becomes 

p(d)  =  E{p(r,t)  p*(r+d,  t)} 

M  M 

=  Z  Z  an  am  exP^[(kn-km)r"kmd)]}  9> 

n=l  m=l 

where  p^n.m)  =  E{exp(  j<j>n)  exp(-j<J>m) } . 

For  uncorrelated  modes  we  get  p^(n,m)  =  0,  n^m,  and  hence 
M 

p(d)  =  E  a*  exp(-jk  d) .  (Eq.  10) 

n=l  n  n 

In  this  case  the  field  is  independent  of  the  absolute  distance  between 
source  and  receiver  r,  except  for  modal  attenuations  contained  in  an>  It 
depends  only  on  the  spacing,  d,  between  sensors.  In  particular,  for  a 
horizontal,  equally-spaced  array,  d  has  to  be  replaced  by 

dik  =  (i  -  k)  cose  dQ.  (Eq.  11) 

In  this  case  the  covariance  matrix  across  the  array  is  Toeplitz.  Except 
for  Sect.  3.5  the  following  investigations  are  based  on  a  signal  model  as 
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given  by  Eqs.  10  &  11.  Additional  white  noise  is  assumed.  The  following 
parameters  are  kept  constant  if  not  mentioned  separately: 


winter  sound-speed 

profi 

ile 

as 

in  Fig. 

la 

summer  sound-speed 

profi 

ile 

as 

in  Fig. 

lb 

water  depth 

100 

m 

bottom- layer  depth 

1 

m 

array  depth 

50 

m 

frequency/spacing 

75 

Hz 

/ 

5  m 

200 

Hz 

/ 

20  m 

900 

Hz 

/ 

50  m 

signal-to-noise  ratio  20  dB 
number  of  sensors  5 


SOUNO  SI'EtD  (m/s) 

1500  15/0  1540  1560  1580  1600 


FIG.  1  SOUND-SPEED  PROFILES 

a)  winter  (SVP1) 

b)  summer  (SVP2) 
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2  METHOD 

Supposing  a  pilot  source  with  known  bearing  transmitting  test  signals,  how 
can  the  signals  received  by  the  array  be  used  to  calculate  a  set  of 
steering  matrices  H(B)  so  as  to  keep  bearing  errors  low?  An  obvious 
solution  is  to  move  the  pilot  source  in  a  circle  and  to  store  test-signal 
covariance  matrices  for  any  bearing  of  interest.  This  involves,  of  course, 
that  for  any  bearing  of  interest,  e.g.  every  \  degree,  the  test  source  has 
to  be  located  by  other  than  acoustic  methods,  which  is  a  tedious  procedure. 

The  purpose  of  the  following  is  to  find  a  method  for  low-error  bearing 
estimation  given  just  one  test  signal  arriving  from  the  known  direction 
The  idea  is  to  approximate  the  actual  covariance  function  achieved  from 
a  test  source  by  a  covariance  function  of  the  following  form: 

pa(i-k>  6t)  =  aa(i-k,  Bt)  Pa(i-k,  Bt),  (Eq.  12) 

where 

aa(i-k,  Bt)  =  sinc[-lp  *u*d0(i-k)cosBt]  (Eq.  13) 

is  the  envelope  and 

Pa(i-k,  Bt)  =  exp[j|^  *v*d0(i-k)cosBt]  (Eq.  14) 


is  the  phase  term.  In  other  words,  the  measured  test  covariance  function 
Pt(i-k,  Bt),  describing  a  certain  spatial  distribution  of  the  test  signal, 
has  to  be  approximated  by  a  complex  function  of  the  kind  of  Eq.  12,  which 
describes  a  spatial  interval  with  uniform  energy  distribution.  Thereby  the 
unknown  energy  distribution  of  the  test  signal,  which  is  characteristic  for 
the  actual  state  of  the  channel,  is  replaced  by  an  equivalent  angle 
interval  of  width  u  (difference  between  two  directional  cosines)  and  the 
mid-direction  v  cosB^)-  The  approximation  has  to  be  carried  out  separately 
for  amplitude  and  phase  term  by  maximizing  the  generalized  cos2(*) 
between  the  vectors  of  the  test  signal  and  those  containing  the  terms  of 
Eqs.  13  and  14: 


al  a  aT  a 
max  — t  — a  — t  —a 


it  * t  h 


P*  P  P*  P 

max  -t  -a  -t  -a 

V  P?  p*  EJ  K 

— t  — l  — a  — a 


(Eq.  15) 


where  a*  and  are  the  vectors  of  amplitudes  and  phase  coefficients  of  the 
measured  test  covariance  matrix  and  a^  and  are  vectors  containing  Eqs. 

13  and  14  respectively,  i.e.  amplituHe  and  pnase  coefficients  of  the 
approximated  covariance  function.  Maximization  of  expressions  15  yields 
the  unknown  parameters  u  and  v  of  the  equivalent  angle  interval  with 
uniform  energy  distribution.  A  simple  algorithm  for  maximizing  expressions 
15  is  shown  in  Fig.  2. 
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The  next  step  is  to  test  how  good  the  approximation  is,  or,  in  other  words, 
how  far  the  assumption  of  uniform  energy  distribution  is  valid.  This  is 
done  by  investigating  the  covariance  values  calculated  by  means  of  Eq.  12 
and  using  the  values  u  and  v  achieved  from  the  approximation  procedure  in  a 
square  matrix  H(8);  furthermore,  H(8)  has  to  be  inserted  in  one  of  the 
generalized  resolution  methods  (Eqs.  2,5,6).  If  there  is  no  difference 
between  the  known  bearing  of  the  test  source  and  the  position  of  the 
maximum  obtained  from  the  processor  we  conclude  that  the  energy  of  the  test 
signal  is  roughly  uniformly  distributed,  i.e.  the  modal  amplitudes  a  in 
Eq.  10  are  equally  spaced  and  constant  for  all  n.  n 

If  the  assumption  of  uniform  energy  distribution  is  not  verified  we  get  a 
certain  bearing  error  by  steering  the  approximated  beamformer  matrix  H(8) 
over  the  test  source;  this  can  be  used  for  correcting  all  bearings 
according  to 

Bc  =  3h  -  <S(cos8h  /  cosBt)a.  (Eq.  16) 

Here  6h  is  the  direction  at  which  the  beamformer  matrix  H  is  looking,  8t  is 
the  bearing  of  the  pilot  source,  8C  is  the  corrected  bearing,  and  6  is  the 
observed  bearing  error.  This  correction  is  based  on  the  idea  that  as  the 
width  of  the  bearing  interval  is  proportional  to  cosB,  the  bearing  error 
will  be  proportional  to  cosB  as  well.  For  instance,  for  the  steering 
matrix  pointing  in  the  direction  of  the  test  source  one  gets  6h  =  3t  and 
hence  8^  =  8h  -  6,  which  is  in  fact  the  appropriate  correction.  The 
exponent  a  is  to  make  the  correction  somewhat  nonlinear  in  order  to  cope 
with  the  nonlinear  distortion  of  the  modal  energy  distribution  for 
different  values  of  8.  For  instance,  a  source  located  at  an  angle  80  that 
causes  a  certain  angle  interval  to  be  uniformly  covered  by  signal  energy 
will  never  appear  to  be  uniformly  distributed  at  any  other  bearing  6+8  ; 
this  is  due  to  the  nonlinear  distortion  of  the  bearing  scaling.  The  ° 
choice  of  a  is  not  very  critical;  a  =  5  to  15  usually  yields  satisfactory 
results.  An  optimum  value  cannot  be  given  because  it  depends  on  the  actual 
amplitude  distribution,  which  depends  again  on  the  actual  channel  parameters. 
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3  NUMERICAL  EXAMPLES 

This  chapter  discusses  several  plots  (Figs.  3  to  14)  that  demonstrate  how 
the  method  outlined  above  works  and  what  the  limitations  are.  Throughout 
the  examples  the  maximum  entropy  method  (Eqs.  4  &  6)  is  used  for  bearing 

estimation.  On  top  of  each  plot  are  two  subfigures  that  indicate  the  true 
★ 

position  ($)  and  the  modal  powers  a^  in  logarithmic  scale.  The  upper 
subfigure  always  denotes  the  pilot  source  used  to  estimate  the  parameters 
u  and  v  (spatial  spread  and  mean  direction).  The  lower  one  represents  the 
target  source  whose  bearing  is  to  be  estimated.  The  lower  one  represents 
.the  target  source  whose  bearing  is  to  be  estimated.  There  are  always  four 
curves.  The  first  ("TEST")  is  achieved  by  steering  the  processor  (Eq.  6) 
containing  the  approximated  steering  matrix  (Eq.  12)  over  the  test  source. 
This  curve  shows  how  well  the  assumption  of  uniform  energy  distribution 
fits  the  actual  channel  considered.  Furthermore,  we  get  the  bearing  error, 
which  is  used  for  compensation.  The  OMEM  (optimum  maximum  entropy  method) 
curve  shows  the  response  of  the  processor  (Eq.  6)  for  H(B),  which  is  the 
exact  set  of  signal  covariance  matrices  as  modelled  by  the  normal -mode 
program,  (e.g.  Eqs.  9  &  10),  and  hence  is  the  optimum  bearing  estimator. 

The  AMEM  curves  (approximated  maximum  entropy  method)  represent  the 
response  of  the  processor  (Eq.  6),  with  H(B)  being  calculated  by  means  of 
Eq.  12  using  the  values  u,  v  and  6  achieved  from  the  approximation 
procedure  of  Fig.  2.  Finally,  the  MEM  curve  shows  the  power  response  of 
the  conventional  maximum  entropy  method  as  given  by  Eq.  4. 


3.1  Different  Channels  and  Frequencies 

Figures  3  to  8  show  examples  for  several  sets  of  parameters:  the  winter 
and  summer  sound-speed  profiles  of  Fig.  1,  frequencies  of  75,  200,  and  900  Hz, 
and  source  depths  of  2,  50,  and  96  m.  The  test  source  is  always  located  at 
10°,  the  target  at  20°. 

The  first  impression  is  that  the  conventional  MEM  yields  high  resolution, 
but,  in  most  cases  gives  considerable  bearing  errors.  Ir.  particular,  the 
MEM  attempts  to  resolve  the  internal  structure  of  the  channel  response  to 
the  source,  which  is  not  very  meaningful  for  source  location.  In  some 
cases  single  modes  are  resolved  (e.g.  Fig.  3b).  The  optimum  method  (OMEM) 
points  always  at  the  true  target  position,  as  expected.  The  method 
suggested  in  this  paper  sometimes  achieves  the  right  bearing  (Figs.  4a, 

8a);  in  all  other  cases  bearing  errors  of  0.5°  to  1.5°  are  observed.  It  is 
noticed  that,  except  for  a  few  examples  (Figs.  3a,  3c,  6a,  7a),  the  bearing 
estimate  is  upward  biased,  i.e.  the  bearing  estimate  is  higher  than  the 
true  bearing  of  the  target.  The  reason  is  that  in  these  examples  the 
bearing  error  is  slightly  over-compensated  by  the  term  of  Eq.  16.  This 
offset  can  be  reduced  by  choosing  the  exponent  a  to  be  of  a  higher  value. 

In  the  examples  of  Figs.  3  to  8  we  have  chosen  a  =  5.  The  choice  a  =  10  to 
15  would  yield  even  better  results  on  the  average  than  those  presented 
here. 
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Apart  from  the  fact  that  in  most  of  the  examples  presented  here  the 
conventional  resolution  method  (MEM)  yields  two  peaks  instead  of  one,  we 
observe  that  there  is  almost  no  example  indicating  that  under  certain 
conditions  the  MEM  is  superior  to  the  AMEM  (compare  the  bearing  of  the 
highest  peak  of  the  MEM  curves  with  the  AMEM  response).  Only  in  Fig.  3b  is 
the  bearing  of  the  higher  MEM  peak  closer  to  the  true  bearing  than  the 
AMEM.  There  are,  however,  several  examples  (e.g.  Figs.  3a,  3c,  6a,  7a,  8a) 
that  show  a  significant  improvement  of  the  bearing  estimate  achieved  by  the 
AMEM.  The  AMEM  is  superior  in  particular  if  the  major  part  of  the  energy 
is  in  the  higher  order  modes,  which  happens  if  the  source  is  close  to  the 
surface  (e.g.  Figs.  3a,  3c,  7a,  8a). 


3.2  Aperture 

Figure  9  uses  the  example  shown  in  Fig.  8a  to  investigate  the  upper  limit 
of  resolution  of  the  AMEM.  Different  plots  show  the  responses  of  arrays 
with  increasing  spacing  (2,  4,  8,  12  m),  corresponding  to  1.19,  2.38,  4.76, 
7.14A.  There  are  always  five  sensors. 

Figure  9a  shows  that  (for  the  example  chosen)  up  to  a  6A  aperture  the  AMEM 
is  identical  to  the  optimum  processor  (OMEM).  The  maximum  resolution 
achievable  is  given  by  using  a  12A  aperture  (Fig.  9b).  A  slight  difference 
between  OMEM  and  AMEM  response  can  be  noticed.  Doubling  the  spacing  once 
more  to  24A  (Fig.  9c)  yields  already  a  considerable  degradation  of 
resolution  and  deviation  from  the  OMEM  response.  The  first  grating  lobes 
can  be  seen  on  the  right. 

Figure  8d  shows  the  performance  of  an  array  with  a  36A  aperture.  It  can  be 
recognized,  in  particular,  that  contrary  to  the  examples  discussed  before, 
the  TEST  curve  (steering  the  uncorrected  processor  over  the  pilot  source) 
is  completely  offset.  Consequently  the  corrected  AMEM  processor  also 
points  in  the  wrong  direction.  The  reason  for  this  behaviour  is  simply 
that  the  approximation  algorithm  (Fig.  2)  did  not  converge  properly,  which 
may  of  course  lead  to  any  erroneous  result.  If  the  approximation  has 
failed,  no  reasonable  result  can  be  expected  from  the  AMEM. 

As  a  conclusion,  we  should  first  keep  in  mind  that  resolution  is  always 
limited  by  any  lack  of  a  priori  knowledge  about  the  signal  to  be  resolved. 
According  to  the  performance  of  the  AMEM  this  means  that  the  resolution 
achieved  depends  on  how  the  measured  covariance  function  can  be  approxi¬ 
mated  by  the  complex  sinc-function  (Eq.  13).  Obviously,  the  approximation 
works  better  for  small  apertures  than  for  larger  ones;  hence  there  is  a 
natural  limit  for  resolution.  Further,  we  can  say  that  a  generalized 
beamformer-type  processor  (Eq.  2)  using  the  approximation  procedure  to 
calculate  H(8)  is  not  at  all  applicable  for  bearing  estimation  because  it 
needs  a  very  large  aperture  in  order  to  achieve  resolution  comparable  to 
the  adaptive  methods  (MLE,  MEM),  which  prevents  convergence  of  the 
approximation  algorithm. 

As  it  turns  out  that  the  approximation  procedure  is  limited  to  rather  small 
apertures  (about  12A)  the  general  policy  should  be  to  use  a  large  array  and 
conventional  or  adaptive  beamforming,  or  any  other  kind  of  detection  scheme 
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[5]  based  on  the  CBF  that  maximizes  the  signal-to-noise  ratio  in  order  to 
detect  targets  from  noisy  background.  As  a  second  step,  only  a  small  part 
of  the  array  should  be  used  for  bearing  estimation  using  the  AMEM. 


3.3  Position  of  Pilot  Source  and  Target 

This  section  looks  at  the  questions  of  whether  there  is  a  preferable  pilot- 
source  position  and  if  one  pilot  source  at  one  known  bearing  is  sufficient 
to  model  the  steering  matrix  H(£)  for  the  whole  bearing  range  of  interest. 

Figure  10  shows  the  same  examples  as  before,  but  with  the  pilot  source 
being  moved  around.  As  can  be  seen,  the  bearing  estimate  is  quite 
independent  of  the  position  of  the  test  source.  Of  course,  for  bearing 
angles  close  to  90°  (broadside)  no  good  estimates  can  be  expected  for  the 
values  u  and  v  in  Eqs.  12,  13  and  14,  because  the  angle  interval  of  signal 
arrivals  approaches  zero.  The  method  seems  also  to  be  very  insensitive  to 
any  change  in  target  position,  as  seen  in  Fig.  11. 

Figure  12  shows  two  examples  in  which  the  pilot  source  is  located  at 
endfire  direction.  The  examples  chosen  are  the  same  as  in  Figs.  7a  and  7c, 
except  that  instead  of  the  pilot  source  being  at  10°  it  is  here  at  endfire. 
It  is  seen  that  in  Fig.  12a  the  AMEM  is  as  good  as  it  was  in  Fig.  7a. 

Figure  12b,  however,  shows  that  there  are  some  situations  in  which  the 
endfire  position  is  not  useful  for  positioning  the  pilot  source.  Fig.  7c 
showing  a  much  better  bearing  estimate.  This  can  be  explained  by  the 
spatial  aliasing  effect  at  endfire.  The  channel  response  to  a  source  close 
to  endfire  is  corrupted  by  its  reflected  image  about  the  value  of  B  =  0°. 
This  aliased  signal  response  is  in  many  cases  useless  for  predicting  the 
spatial  spread  and  mean  direction  of  a  source  at  angles  different  from 
endfire  and,  hence,  much  less  aliased. 


3.4  Mismatch  in  Depth 

So  far  it  has  been  assumed  that  the  depths  of  the  test-source  and  the 
target  are  the  same.  In  the  two  examples  in  Fig.  13a  &  b  it  is  assumed 
that  the  target  is  at  2m  depth  and  that  the  pilot  source  is  at  30  m  and 
50  m  depths  respectively.  Except  for  the  changes  in  pilot  source  depth, 
these  examples  are  the  same  as  that  given  in  Fig.  7a.  These  examples  show 
that  there  can  be  considerable  differences  in  energy  distribution  between 
test  signal  and  target.  In  some  instances  it  may  happen  that  the  bearing 
estimate  is  corrupted  seriously  by  mismatch  between  pilot-source  response 
and  response  to  the  target,  as  shown  in  Fig.  13a.  In  other  instances,  such 
as  in  Fig.  13b,  a  good  estimate  is  achieved  although  the  spatial  character¬ 
istics  are  quite  different.  A  general  statement  about  the  sensitivity  to 
mismatch  in  depth  between  pilot  source  and  target  cannot  be  made. 

Therefore,  if  possible,  the  pilot  source  should  be  kept  at  the  target  depth 
of  interest. 
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3.5  Correlation  among  Modes 


So  far  it  was  assumed  that  modes  are  uncorrelated  among  each  other. 


=  1 
=  0 


n  =  m 
n  ^  m. 


i  .e. 


In  this  case  wavefronts  are  random,  non-planar,  and  range  independent  in 
the  mean.  The  spatial  covariance  matrix  is  Toeplitz  and,  hence,  the  best 
match  between  the  approximation  method  and  the  field  is  achieved  because 
the  steering  matrix  H(8),  as  calculated  by  Eqs.  12,  13  and  14,  is  Toeplitz. 

If,  however,  there  is  some  correlation  between  modes  the  field  becomes 
range  independent  and  the  covariance  matrix  non-Toeplitz.  If,  for  instance, 

3<J>n 

E{e  e  }  =  1  for  all  n,m  one  receives  non-random  non-planar 
wavefronts  whose  shapes  depend  on  the  distance  between  source  and  receiver. 
Therefore  it  can  be  expected  that  the  resolution  of  any  processor  (except 
for  the  OMEM,  which  contains  a  priori  knowledge  about  the  modal  amplitudes 
and  phases)  is  degraded  due  to  mismatch  to  the  field  (mismatch  of  the 
steering  vector  as  suboptimum  steering  matrix  to  non-planar  wavefronts). 


In  other  words,  if  the  modes  are  uncorrelated,  i.e.  phase  among  modes  are 
entirely  random,  then  the  signal  energy  is  distributed  in  the  mean  as  shown 
in  the  subfigures  of  the  plots.  If  modes  are  correlated,  the  shapes  of  the 
wavefronts  are  distorted  additionally  by  mode  interference,  which  may  cause 
a  significant  deviation  from  a  rectangular  energy  distribution  that  the 
AMEM  is  based  on.  As  a  consequence,  the  signal  field  cannot  be  approximated 
sufficiently  well  by  the  sinc-function  (Eqs.  12,  13  &  14),  thus  causing 
additional  mismatch  between  the  AMEM  processor  and  the  signal  field. 


We  use  the  following  correlation  model 


p4>( n »m)  =  P0(n’m)  P0(n,m)  >  0 

=  0  pQ(n,m)  <  0,  (Eq.  17) 

where 

/  m  m  \  1  C I  n-m  I 

P0(n»m)  =  1  -  'm  -1  • 

For  C  =  0  we  get  p0(n,m)  =  1,  i.e.  modes  are  entirely  correlated, 

C  >  M  yields  uncorrelated  modes.  Some  results  are  shown  in  Fig.  14.  It 
can  be  seen  that  for  C  =  1  (Fig.  14a)  (triangular  correlation  among  modes) 
rather  high  resolution  is  achieved.  The  AMEM  obviously  gives  a  rather  good 
bearing  estimate  (0.5°  offset).  The  correlation  is  still  weak  enough  for 
appropriate  approximation  of  the  signal  covariance  matrix  by  the  sinc- 
function  of  Eq.  12.  For  C  =  0.1  (Fig.  14b)  the  correlation  function  is  a 
triangle  on  a  pedestal  of  0.9  height.  A  considerable  loss  in  resolution  is 
observed  and  the  mismatch  of  the  steering  matrix  (in  the  case  of  the  AMEM) 
or  the  conventional  beamformer  (in  the  case  of  MEM)  now  becomes  significant. 
With  total  correlation  (C  =  0,  Fig.  14c)  the  resolution  achieved  is  even 
worse  (except  for  the  OMEM,  which  is  perfectly  matched  in  any  case). 

However,  in  all  cases  considered,  the  AMEM  appears  to  be  superior  to  the 
conventional  MEM. 
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5  EXAMPLES  FOR  75  Hz  AND  WINTER  SOUND-SPEED  PROFILE 

a)  2  m  pilot-source  depth 

b)  50  m  pilot-source  depth 

c)  96  m  pilot-source  depth 
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FIG. 


EXAMPLES  FOR  900  HZ  AND  WINTER  SOUND-SPEED  PROFILE 

a)  2  m  pilot-source  depth 

b)  50  m  pilot-source  depth 

c)  96  m  pilot-source  depth 
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FIG.  8  EXAMPLES  FOR  900  Hz  AND  SUMMER  SOUND-SPEED  PROFILE 

a)  2  m  pilot-source  depth 

b)  50  m  pilot-source  depth 

c)  96  m  pilot-source  depth 
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11  EXAMPLE  OF  FIG.  8  BUT  WITH  CHANGES  IN  THE  BEARING 
OF  THE  TARGET 

a)  2  m  pilot-source  depth 

b)  50  m  pilot-source  depth 

c)  96  m  pilot-source  depth 


OUTPUT  LEVEL  i OB 5  OUTPUT  LEVFL  l OB 

,-4Q . 00  -30.00  _ -20  oo  JO.  00  0^.00 _  10-00  20-00  o-4Q.Q0  -30 ■  00  -20-00  -lQ.OO  0^ 


SACLANTCEN  SR-40 


S-00 


10-00 


T^ToO  20  do  2T.Q0  30 - 00  3^-00  40.00 

BEARING  COEGJ 


45.00 


LES  WITH  THE  PILOT  SOURCE  AT  ENDFIRE  (200  Hz,  SUMMER  PROFILE) 

2  m  pilot-source  depth 
96  m  pilot-source  depth 


SACLANTCEN  SR-40 


FIG.  13  EXAMPLES  WHERE  PILOT  SOURCE  AND  TARGET  ARE  NOT  AT 
THE  SAME  DEPTH 

a)  source  at  2  m,  target  at  30  m 

b)  source  at  2  m,  target  at  50  m 
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CONCLUSIONS 


A  method  (AMEM)  has  been  described  to  reduce  the  bearing  error  of  hori¬ 
zontal  line  arrays  caused  by  multipath  effects  of  sound  propagating  in 
shallow  water.  The  procedure  consists  of  the  following  steps: 

1.  Use  a  pilot  source  of  a  certain  known  bearing  angle  between  10° 
and  45°  (0°  =  endfire)  or  135°  and  170°  and  estimate  the  spatial  correla¬ 
tion  coefficients  for  all  delays  of  the  array,  i.e.  p^-j,  p-|2,  *  *  ,plN' 

2.  Approximate  the  measured  correlation  function  by  a  sinc(*)exp(j(*))- 
function  (Eqs.  12,  13  &  14)  using  the  algorithm  of  Fig.  2  and  achieve  the 
values  u  (width  of  equivalent  angle  interval)  and  v  (mid  direction). 

3.  Calculate  steering  matrices  H(B)  by  means  of  Eqs.  12,  13  &  14 
using  the  estimated  values  of  u  and  v. 

4.  Insert  H(8)  in  a  high  resolution  method  (Eqs.  5  &  6),  preferably, 
steer  the  processor  over  the  pilot  source,  and  determine  what  bearing  error 
is  characteristic  of  the  existing  channel  conditions. 

b.  Use  the  bearing  error  6  for  bearing  correction  of  H(B)  by  means 
of  Eq.  16. 

6.  If  there  is  some  directional  noise  present  when  the  steering 
matrix  is  being  updated,  then 

a.  estimate  spatial  noise  covariance  function  n,,  .n-jp- •  •niw, 
when  no  test  signal  is  present, 

b.  estimate  the  spatial  covariance  function  r,, »ri2”  • •riw» 

when  signal  and  noise  are  present,  c 

c.  subtract  the  noise  component  from  the  measured  covariance 
function  Pll  =  r^-n^,  p]2  =  r]2-n12,  p]N  =  r]N-n1fr 

The  method  has  been  investigated  numerically  by  application  to  a  normal¬ 
mode  sound-propagation  model,  assuming  that  modes  are  uncorrelated.  The 
results  can  be  summarized  as  follows: 

a.  The  AMEM  response  to  a  point  source  is  usually  one  peak  in 
bearing.  Conventional  resolution  methods  like  MEM  or  MLM  may  yield  a 
thousand  or  more  peaks,  depending  on  the  spatial  energy  distribution  of  the 
signal  and  the  aperture  of  the  array. 

b.  The  method  has  been  tested  for  different  conditions  of  the 
channel  (two  sound-speed  profiles,  three  frequencies,  three  source  depths). 

In  most  cases  observed,  the  maximum  of  the  AMEM-response  was  closer  to  the 
true  bearing  of  the  target  than  the  maximum  peak  determined  by  the  MEM. 

c.  The  maximum  aperture  to  be  used  is  about  12A.  Beyond  12A  the 
algorithm  of  Fig.  2  does  not  converge  sufficiently  well,  thus  causing  a 
mismatch  between  steering  matrix  and  field  that  leads  to  degradation  of 
resolution. 
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d.  The  method  is  useful  for  bearing-error  reduction  for  targets 
located  between  10°  and  170°.  Targets  at  about  90°  (broadside)  do  not  of 
course  need  to  be  corrected.  At  endfire  direction  (0°  to  10°)  the  method 
is,  in  many  cases,  corrupted  by  spatial  aliasing  and  should  be  used 
carefully. 


e.  The  pilot  source  should  be  located  at  angles  between  10°  and  45°. 
The  accuracy  of  the  AMEM  depends,  of  course,  on  how  precisely  the  position 
of  the  pilot  is  known. 

f.  The  AMEM  (or  AMLM)  is  superior  to  conventional  adaptive  methods 
(MEM,  MLM) ,  particularly  if  the  major  part  of  the  energy  is  contained  in 
the  higher-order  modes,  i.e.  for  a  source  located  close  to  the  surface  or 
to  the  bottom. 

g.  Even  if  the  field  is  inhomogeneous  due  to  mode  interference  the 
AMEM  seems  to  be  superior  to  the  conventional  methods. 

h.  The  si nc-f unction  seems  to  be  an  appropriate  way  of  approximating 
the  envelope  of  the  covariance  function  because  the  assumption  of  uniform 
energy  distribution  is  available.  Nevertheless,  in  certain  cases,  other 
functions  might  be  even  better,  which  can  always  be  judged  by  consideration 
of  the  approximation  error  1-e  in  the  algorithm  of  Fig.  2.  Of  course, 
whatever  function  is  used,  it  has  to  be  one  that  depends  only  on  u*cosB  in 
order  to  adjust  it  to  the  interval  width  varying  with  bearing. 

i.  The  AMEM  is  basically  sensitive  to  mismatch  between  pilot-source 
depth  and  target  depth.  Therefore  the  pilot  source  should  be  located  at 
the  depth  of  interest  for  target  location. 

j.  The  method  outlined  above  is  adaptive  in  that  it  allows  the 
steering  matrix  to  be  updated  with  changing  environmental  conditions. 

k.  Basically,  the  method  is  based  on  the  assumption  that  there  are 
phase  fluctuations  among  modes,  which  is  supposed  to  be  valid  for  higher 
frequencies  rather  than  for  lower  ones.  The  validity  of  this  assumption 
has  to  be  proved  experimentally. 
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